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Abstract 

Discrete symmetry breaking and possible restoration at finite temperature T 
are analysed in 2D Gross-Neveu model by the real-time thermal field theory in the 
fermion bubble approximation. The dynamical fermion mass m is proven to be scale- 
independent and this fact indicates the equivalence between the fermion bubble dia- 
gram approximation and the mean field approximation used in the auxialiary scalar 
field approach. Reproducing of the non-zero critical temperature Tc = 0.567m(0), 
(m(0) is the dynamical fermion mass at T = 0), shows the equivalence between 
the real-time and the imaginary-time thermal field theory in this problem. How- 
ever, in the real-time formalism, more results including absence of scalar bound 
state, the equation of criticality curve of chemical potential-temperature and the 
ln(rc/r) behavior of at T ~ can be easily obtained. The last one indicates 
the second-order phase transition feature of the symmetry restoration. 
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1 Introduction 

In 4D model of NJL-form [1] and 3D Gross-Neveu (GN) model [2], by means of analyses 
of dynamical fermion mass in the real-time thermal field theory in the fermion bubble 
approximation, we have acquired a quite clear understanding of restorations of the spon- 
taneously broken symmetries at high temperature and their second-order phase transition 
feature [3,4]. A natural idea is to extend the same analysis to 2D GN model. However, 
once entering one-space dimension, one will find some special features. By the Mermin- 
Wagner- Coleman Theorem [5], a continuous symmetry can never be spontaneously broken 
in 2 dimensions. Consequently we will have to confine ourselves to possible spontaneous 
breaking of only some discrete symmetries. The 2D GN model is perturbatively renor- 
malizable and asymptotically free and this is apparently different from the perturbatively 
non-renormalizable 4D NJL model and 3D GN model. I3ased on the auxialiary scalar field 
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approach, many researches on the behavior of 2D GN model at finite temperature have 
been made [6,7,8], though all of them in the imaginary-time formalism of thermal field 
theory. Considering the above facts, in this paper, we will analysed 2D GN model in the 
real-time thermal field theory in the fermion bubble approximation. Through comparison 
between the results obtained by our approach and the known ones, we will be able to 
check the equivalence of the real-time and the imaginary-time formalism of thermal field 
theory in this problem, find the advantages of the real-time formalism and the limitations 
of the fermion bubble approximation and , in the meantime, get deeper insight into the 
2D GN model itself as well. 

The paper is arranged as follows. In Sect. 2 we will analyse the symmetries of 2D GN 
model, mainly its discrete symmetries and their possible spontaneous breaking induced by 
dynamical fermion mass. In Sect. 3 the behavior of the model at temperature T = will 
be outlined. This includes derivation of the renormalized gap equation and the running 
coupling constant in dimension regularization scheme and the proof of that scalar bound 
state could not exist in this model. In Sect. 4 we extend the gap equation to the case of 
T 7^ 0, determine the critical temperature of symmetry restoration and derive the critical 
behavior of the dynamical fermion mass. Finally in Sect. 5 we indicate the advantages 
and the limitations of our approach through the comparison with the known results and 
reach corresponding conclusions. 



2 Symmetry analysis of the model 

The Lagrangian of 2D GN model is expressed by 

N _ o ^ - 

L = J2i>Hrd,^k + f E[^Vfe]'- (2.1) 

k=l ^ fc=l 

where ip = ipk{t,x) are the spinor fields with "color" components. In 2 dimesions, 
the four-fermion coupling constant g is dimensionless thus the theory is perturbatively 
renormalizable. The 7^(/i = 0, 1) and 75 = T^T"*^ are 2x2 matrices which are related to 
the Pauli matrices (7*(i = 1, 2, 3) by 



i.e. we take the "Weyl" representation. In this representation, il}k will be two-component 
complex spinors. We will concern only the discrete symmetries of the model, since con- 
tinuous symmetries can not be spontaneously broken in 2 dimensions. It is easy to see 
that the Lagrangian (2.1) is invariant under the discrete transformation i?, 

V^(i,x)-^-V'(i,a;) (2.3) 
and the ordinary parity V and the time reversal T, 

V^(t, a:) ^7X^,-0:), 

V'(t,x)^7V(-^,^)- (2.4) 
In addition, it is also invariant under the discrete chiral transformation xdi 

ilj{t,x)^^5^{t,x) (2.5) 

and the special parity Vi, 

^{t,x) ^-f^'ip{t,-x) (2.6) 
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We indicate that the fermion mass term —mifjilj will keep R, V and T invariant but not 
Xd and Vi invariant. Therefore, the dynamical generation of the fermion mass will imply 
spontaneous breaking of the discrete chiral symmetry xd and the special parity Vi. The 
spontaneous breaking of xd and Vi and their possilDle restoration at finite temperature 
T is just the problem we will discussed. 



3 The behavior of the model at T = 



The renormalization operation will be conducted in the dimension regularization approach 
[9]. For this purpose, we change the time-space dimension from 2 into 



D = 2-2e 



and the Lagrangian C in Eq.(2.1) into 



N _ N _ 

fc=l ^ k=l 



(3.1) 



(3.2) 



where M is a scale parameter with mass dimension, the mass dimensions [g\ = 0, 
[■0jt] = [■0*^] = (-D — l)/2 and the dimension of the representation of the matrices 
becomes 2^^^. The Feynman rule of the four-fermion vertex will be igM"^'^ . 

Assuming that the four-fermion interactions will lead to the condensates J2k=i{'^'''^k) 
0, then the non-zero dynamical mass 



N 



m(0) = -^M^-^^(V;V,) 



(3.3) 



k=l 



will be generated and the discrete symmetries Xd and Vi will be spontaneously broken. 
After the Wick rotation and introduction of the denotations Z° — — /*, we obtain 

from Eq.(3.3) the gap equation 



1 = ^A^2^/2m2-^|_ 
= gN2''/^M'' 



(27r)^/2 + m2(0) 
_^ r(l-D/2) 
(47r)^/2[m2(0)]i-^/2' 



where we have used the integral formula [10] 

d^l 1 



i-oo (2^' 



V{A-D/2) 



(27r)^ [[2 + m2(0)]^ (47r)^/2r(A)[m2(0)]^-D/2' 
Considering Eq. (3.1), the gap equation can be reduced to 



27r 



1 , 27rM2 ^ 

--^ + \n—— + 0{e) 
e m^{0) 



(3.4) 



(3.5) 



(3.6) 



where 7 = 0.5772 is the Euler constant. The divergent 1/e term can be removed by 
introduction of the counter-term of the four-fermion interactions 



N 



He,. = f M^-^(Z4 - 1) ElV^'V'.)'- 
^ fe=l 



(3.7) 
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The renormalization constant Z4 is defined by the minimal subtraction scheme [9], i.e., 
to one-loop order, 

Z. = l-fi, (3.8) 

then the renormalized gap equation will take the form 

gN 



1 = 
27r 



(3.9) 



Eq. (3.9) gives that 

m^(0) = 27rM^ ^'^^^~~gN ~ ^^'"^^^ 

Since 2D GN model is an asymptotically free theory, as will be argued immediately, the 
g in Eq. (3.10) should be replaced by corresponding running coupling constant. The 
running coupling can be derived by the general renormalization group method. In fact, 
after introducing the counter-term £c.t., we can define the bare four-fermion coupling 
constant 

go = gM^-^'Z^ (3.11) 
which should be independent of the renormalization scale parameter M, i.e. 

m = (3.12) 

Considering Eq. (3.8), Eq. (3.12) will lead to 

dg N ^ , , 

M-^ = / 3.13 

and furthermore, 

1 AT M 1 , , 

— — - = — In h— r 3.14 

g{M) TT Mo g{Mo) ^ ' 

Eq. (3.14) indicates that the running couphng g{M) when M — > 00, i.e. the theory 
is ultraviolet stable and asymptotically free. Now we replace the coupling constant g in 
Eq. (3.10) by the running coupling g{M), then by means of Eq. (3.14), can prove that 

™^(0) = 2.M= exp[-— - ,] = exp[-^^ - tI (3.15) 

i.e. the dynamical fermion mass m(0) is scale-indepcdent, we will obtain the same m(0) 
in the scales M and Mq. On the other words, the order parameter m(0) of the discrete 
symmetry breaking will be independent of the size of the system. 

We can also calculate the renormalized fermionic four-point function T^ip) in scalar 
channel by the methods used in Ref. [11]. The result is 

rf (p) = -Vb' - 4m2(0) + ie\Ko{p), (3.16) 

^ ^ N dx N 1 

Ao(p) = 1— / — 7rr\ ;r~? r = , arctan— ^ ^=, 

477^0 w?{Q) - p^x{l - x) 7rpy4m2(0)/p2 - 1 ^Aw?{Q)/p^-l 

(3.17) 

where we have made the assumption that < < 4772^(0). To examine if = 4771^(0) 
could be a simple pole of T§{p), we may calculate the residual of Ff (p) at p^ = 4m^(0) 
and obtain 

4m2(0) 



lim 



2 -1 



(p) = (3.18) 
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which imphes that (p) does not represent a propogator for any scalar bound state with 
the mass 2m(0). Hence, the scalar bound state does not exist at T = in this model. 
In fact, this conclusion can be generalized to the case of T 0, as will be showed in the 
next Section. 



4 Symmetry restoration and phase transition at fi- 
nite temperature 

We will take the real-time formalism of thermal field theory [12] to research the tempera- 
ture behavior of the dynamical fermion mass. This means that, corresponding to the vac- 
uum expectation value J2k{i^''i^k) replaced by the thermal expectation value J2k {'4^^'4^k)Ti 
we should make the substitution of the fermion propagator 



y— m + is y— m + is 



2TrS{l^ - m')(^+ m) sin^ 0{f, //), 



^""'^^^'^^^ = exp[/3(/o -//)] + ! + exp[/3Ho + ;,)] + !' ^^-^^ 

where ^ is chemical potential and f3 = 1/T. Let m = m{T,^) is the dynamical fermion 
mass at finite temperature T and finite chemical potential /x, then similar to the derivation 
at T = 0, we may obtain the gap equation at T 7^ 



1 _ ^ + In 'tMl -2 [dH S{f - m^) sin^ e{f, pt) 
s J 



(4.2) 



Notice that the integral related to T is convergent and the UV divergence appears only 
in the zero temperature sector. Therefore, by the definition (3.8) of Z4 at T = 0, we 
can remove the 1/s divergence to the one-loop order and obtain the renormalized gap 
equation at T 7^ 



1 = 1^ (in ^ - 7 - 2[/i(y, -r) + h{y, r)] \ (4.3) 



where 



Ii{y,Tr) = J^ ^j^r^ exp( V^^q^ T r) + 1 ^^"^^ 

with the denotations y — m/T and r = /i/T. By the same argument as the one leading 
to Eq. (3.15), it can be proven that the dynamical fermion mass m at finite temperature 
and chemical potential is also scale-independent. This implies that in the bubble fermion 
diagram approximation, the order parameter is always size-independent. 
Substitute the gap equation (3.9) at T = into Eq. (4.3), we will obtain 

^^m(0) ^ (4.5) 

m 

We may assume that there is a critical temperature Tc 7^ so that when T ^ Tc, m ^ 0. 
It is easy to sec that this assumption is consistent with Eq. (4.5). In fact, when m ^ 0, 
in the left-handed side of Eq. (4.5) there will appear In 00 divergence and in the right- 
handed side, the integral at a; — > will approach ^ and has the same In 00 behavior. 
Thus the assumption 7^ is plausible, at least in the bubble diagram approximation. 
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Since when T ~ Tc 7^ 0, m ~ 0, i.e. 7/ -C 1, we can use the high temperature expansion 
of /i(y,^r) [3] and obtain from Eq. (4.5) 



, m(0) 



+ 



Ttt 



-31C(5) + 



7C(3)-93C(5)(|^)2 



1905 



c(7)(^)^; 



(4.6) 



where C(3) = 1.202, ^(5) = 1.037, ... and r/27r < 1 is also assumed. We find that the 
Inm terms in both sides of Eq. (4.6) have cancelled each other. Taking y = i.e. m = 0, 
we will obtain the equation of Tc and the critical chemical potential /^c 



, m(0) 



31C(5)(^)^ 



127C(7)( 



2n' 



+ 



(4.7) 



with the denotation Vr 



Hc/Tc- When chemical potential is zero, Eq. (4.7) gives 

= (e^/7r)m(0) = 0.567 m(0) (4.8) 



This result is identical to the one obtained in the mean field approximation of the auxi- 
aliary scalar field approach in the imaginary-time formalism of thermal field theory [6-8] . 
When the chemical potential is not equal to zero, Eq. (4.7) will be able to determine a 
— Tc criticality curve at rc/2'K < 1. This shows an advantage of the real-time thermal 
field theory approach. Besides this, we may further give the critical behavior of the dy- 
namical fermion mass m i.e. the order parameter of symmetry breaking at T ~ Tc- Since 

1/ -C 1 when T Tc, it is reasonable to keep only the terms of the order C(y^) in Eq. 
(4.6). Then, considering the equation (4.7) of fic — Tc curve (in which fj,c is redenoted by 
11), we may derive the expression for the squared dynamical fermion mass 



m 



T 

ln^ + 
T 



+127C(7)(A) 




/ 



27ry 
7 



- 31C(5) 




C(3) - 93C(5) 



27rJ 





1-1 


V27ry 


/ 7^4 



when T 

If the chemical potential is zero, then Eq. (4.9) gives 



^Tc 



7C(3) T' 



when T^Tc 



(4.9) 



(4.10) 



The logrithmic behavior of for Tc can be compared with the temperature behaviors of 
at T ~ Tc in 4D NJL model and 3D GN model. The latters are [3,4] 



TT 



rp2 _ rp2 



m 



for 4D NJL model 



m 



3 ln(A/r7r) + 7- 1/2' 
' = 8(ln2)r(Tc - T), for 3D GN model 



(4.11) 
(4.12) 



In both cases we have m ~ (Tg — T)^/^ at T ~ Tc and it means that the relevant symmetry 
restorations are the second-order phase transitions. In the case of 2D GN model, based 
on a similar argument to the one given in Ref. [3], we can reach the same conclusion, i.e. 
the discrete symmetry restoration at T > Tc is also the second-order phase transition. 
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Based the method used in Ref. [13], we can obtain the renormalized fermionic four- 
point function Tgrp^p) in scalar channel at finite temperature T which is expressed by 

TItIp) = - + i^) {Kip) + Hip) - iS{p) - R^{p)/[K{p) + H{p) + iSip)]} , 

(4.13) 

where K{p) has the same form as Kq{p) in Eq. (3.17), except that m(0) must be replaced 
by m, and 



p^ -21-p ^ p'^ + 21-p 



(p2 - 2Z • p)2 + £2 (p2 + 2Z • p)2 + £2 



5(/2-m2)sin2^(/°,/x), (4.14) 



AT r 

S{p) = cPl5{f-my[{l+p)^-m^][sm^e{f,fi) cos^ e{f+p°, fi)+cosH{f, fi) sin^ ^(/"V, A*)], 

(4.15) 

N r 

R{p) ^ -jJdH S{P -m^)S[{l+p)^ -m'']sm29{f,p,)sm29{f+p^,fj,). (4.16) 

The form (4.13) of rfy(p) ensures mutual cancellation of the pinch singularities contained 
in Sip) and Rip). It is easy to confirm that when p"^ m^, Hip) is finite but Sip) and 
Rip) are sigular. The latter comes from the feature of the product of the two 5-functions. 
Since Kip) is divergent when p^ ^ rn?, wc will still have the residual of rfj,(p) at p^ = Aw? 

Hence, the conclusion that there is no scalar bound state in the model will be maintained 

at finite temperature. 

Our results about finite temperature behavior of the model are obtained by the real- 
time thermal field theory in the fermion bubble diagram approximation. It should be 
pointed that the fermion bubble diagram approximation has already been able to rcficct 
physical essentiality in 4D and 3D case. However, it has great limitations in 2D case. As 
illustrated in Ref. [8], where the auxialiary scalar field approach was used, for an infinite 
one-space system and a finite iV, in fact there is not symmetry restoration at T 7^ owing 
to the special kink effects when the minimal configuration a^x) of the effective thermo- 
dynamic potential r((Tc(a;), /3) could be a function of space. Only if adx) ~ constant 
is assumed, i.e. in so called the mean field approximation, one could obtain symmetry 
restoration at T 7^ and this is just the result obtained in Refs. [6,7] and in this paper. 
It may be indicated that the scale-independece, or the size-independence of the order 
parameter m(T, /i) proved in this paper is just consistent with the assumption (Jcix) = 
constant in the mean field approximation. Hence we can say that the fermion bubble 
approximation is indeed the mean field approximation named in Ref. [8]. Since the mean 
field approximation is valid only if — > 00 or if iV is large and one consider only a finite 
segment of one-space dimension system, the results about finite temperature obtained in 
this paper are applicable only in the above conditions. 



5 Conclusions 

We have analysed spontaneous breaking and restoration of the discrete chiral symme- 
try xd and the special parity Vi induced by dynamical fermion mass in renormalizable 
and asymptotically free 2D Gross-Neveu model by means of the real-time thermal field 
theory in the fermion bubble diagram approximation. We have proven that the dy- 
namical fermion mass, both m(0) at zero-temperature and m at finite temperature, are 
scale- independent, consistent with the fundamental assumption of the mean field ap- 
proximation in the auxialiary scalar field approach. Hence, the fermion bubble diagram 
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approximation is equivalent to the mean field approximation and must submit to the same 
applicable conditions. All the results, especially reproducing of the non-zero symmetry 
restoration temperature Tc = 0.567 m(0), show the equivalence between the real-time 
and the imaginay-time formalism of thermal field theory in this problem. However, in 
the real-time formalism, one can further prove that there is no scalar bound state in the 
model, easily obtain the equation of criticality curve of chemical potential and temper- 
ature and the ln(Tc/T) behavior of the dynamical squared fermion mass at T ~ Tc. 
The last one clearly indicates the second-order phase transition feature of the symmetry 
restorations. 
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